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^ . Abstract 

We present an exact time-dependent solution to the effective D-brane world- 
,— J . volume theory which describes an inhomogeneous decay of a brane-antibrane sys- 

^ . tem. We compute the quantum energy flux induced by the particle creation in this 

O ■ inhomogeneous and time-dependent background. We find that this calculation is 

essentially equivalent to that of the moving mirror system. In the initial stage, the 
^ . energy flux turns out to be thermal with the temperature given by the inverse of 

j^ ■ the distance between the brane and the antibrane. Later it changes its sign and 

becomes a negative energy flux. Our result may be relevant for the (p)reheating 
process or/and the evolution of cosmic string network after stringy brane inflation. 
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1 Introduction 

The brane-antibrane system [T] has been studied intensively until recently as a typical 
example of unstable non-supersymmetric systems in string theory j2] • When a brane and 
an antibrane are close to each other, the open string spectrum includes tachyons. The 
condensation of open string tachyons offers us an important time-dependent background. 

When the distance between a brane and an antibrane is larger than the string scale, 
there exits no tachyon. However, by an instanton effect or a small dynamical perturbation 
induced by collisions with other objects, it is possible that the distance between them 
becomes the string scale at a particular region or a point. In this setup, we expect that 
there exists a localized open string tachyon and the tachyon condensation occurs only 
in that region. This leads to a recombination of the brane and antibrane. Once this 
happens, this system annihilates via the time evolution and eventually the whole system 
will disappear radiating closed strings. However, the study in this direction has been 
rather few even at present. 

The brane-antibrane system plays important roles not only in theoretical foundations 
of string theory but also in its cosmological applications. (See ref. [3] and references therein 
for a recent review of string cosmology.) Indeed, a string theoretic realization of hybrid 
inflation, called warped brane inflation [1], was made possible by non-trivial dynamics 
of the brane-antibrane system. In this scenario, the distance between a brane and an 
antibrane embedded in a warped geometry plays the role of inflaton. Our 4-dimensional 
universe parallel to the world-volume of the brane-antibrane system expands with an 
accelerated expansion rate while the interbrane distance slowly decreases [5]. When the 
interbrane distance becomes as short as the local string length, a tachyon appears and 
starts rolling. This process of brane-antibrane annihilation releases sufficient energy to 
reheat our universe and starts the hot big-bang cosmology. 

In hybrid inflation in fleld theory [6l [7] , the reheating process begins with the so called 
tachyonic preheating [8]. This highly non-linear, inhomogeneous dynamical process occurs 
due to the tachyonic instability near the top of the water-fall potential. The tachyonic 
instability converts most of the energy into that of colliding classical waves very rapidly, 
within a single oscillation. The tachyonic preheating is a typical example of processes in 
which inhomogeneities are crucial. 

As in the tachyonic preheating of the fleld-theoretic hybrid inflation, inhomogeneities 
are expected to play essential roles also in the reheating process of the warped brane 
inflation. In view of this, the brane-antibrane annihilation should begin locally in re- 
gions where the interbrane distance reaches the string length earlier than the other parts. 



and those regions should turn into "expanding holes" on the world- volume of the brane- 
antibrane system. The whole annihilation process proceeds as those "expanding holes" 
collide with each other and percolate. Therefore, the study of inhomogeneous annihilation 
of brane-antibrane pair is an important subject towards our better understanding of the 
end of the brane inflation and the beginning of the thermal history of our universe. 

Motivated by these, in this paper we would like to analyze the real time evolution of 
a brane-antibrane system after the local annihilation or recombination occurs (refer to 
[9] for an analysis of the process of the recombination). We first present a simple exact 
solution to the DBI action of D-branes which describes the time evolution of a brane- 
antibrane system. Then we will calculate the quantum energy flux produced by the open 
string creation. In particular, we will mainly focus on the brane-antibrane system of a 
D-string and an anti D-string (i.e. Dl — Dl system) for simphcity by a technical reason 
that we can apply the two dimensional conformal field theory. The argument of general 
Dp-branes can also be possible in principle and will not be essentially different. As we will 
show, our computation of energy flux has a close relation to the analysis of the moving 
mirror [TOl [11]. Finally, we will interestingly notice that we encounter a negative energy 
flux in this process of decaying branes. 

The inhomogeneous annihilation of D-string and anti D-string considered in this paper 
may be regarded as a lower-dimensional proto-type of the "expanding hole" mentioned 
above [12]. The process considered in this paper is probably relevant also for the dy- 
namics of cosmic superstring network formed after the warped brane inflation. Indeed, 
the recombination of strings is one of the most important factors for the evolution of the 
cosmic superstring network. (See ref. [13] and references therein for a review of cosmic 
superstrings.) 

This paper is organized as follows. In section two, we present an exact solution which 
describes the decay of brane-antibrane systems. In section three, we study the induced 
metric on this time- dependent world-volume of the decaying brane. In section four, we 
compute the energy flux in this process induced by the particle creation mechanism, 
relating this calculation to the one in the moving mirror system. We also discuss its 
string theoretic interpretations. In section flve we summarize conclusions and discuss 
future problems. 

2 Real Time Decay of D-string anti D-string Pair 

Consider a parallel D-string anti D-string pair separated by a distance d. When 
d > 2txIs in bosonic string (or d > \/2itIs in superstring), the open string between the 



D-string and the anti D-string does not have any tachyonic modes {Ig is the string length). 
Therefore in this case, the decay will be induced by their relative motion. 

Even though we started with the system of a brane and an antibrane, we can regard 
the total system as a single D-string after the recombination occurs. Then we can describe 
the motion of the system by the DBl action of D-string as long as the stringy corrections 
are negligible. 

2.1 Exact Solution of Time Evolution 

We assume that the D-string is included within the (2-1-1) dimensional hypersurface i?^'^, 
specified by the coordinate {t, x, y) in the D = 10 (or D = 26) dimensional superstring (or 
bosonic string) theory. The profile of the D-string is described by the function y = y{t, x). 
Then the Lagrangian of the DBI action for the D-string is written as follows 

L=^l-y^ + y'\ (2.1) 

where y = -^ and y' = -^. 

Let us find its solution by assuming the form y{t,x) = f{x) + g{t). The equation of 
motion for (12.11) leads to 

m+n+f"{9'-^) = o. (2.2) 

Thus we find -r^ = tttj2 =const. This can be solved as follows 

t X 

g(t) = a log cosh— , f{x) = —a log cos— . (2.3) 

a a 

In this way we find the time-dependent solution 

e^/" cos(x/a) = cosh(t/a), (2.4) 

where x runs the range —fa < x < ^a. Actually, its Euclidean version e'^^" cos{x/a) = 
cos{z/a) is known as the Scherk minimal surface in differential geometry. 

This solution (12.41) represents a decay of a D-string and an anti D-string. Indeed, 
when y is very large at a fixed time t, the coordinate x approaches to either of two 
values X = ±|a. These two branches correspond to the D-string and the anti D-string 
separated by the distance d = ira. Since we assume that there are no tachyons between 
them, we require a ^ Ig- Notice also that y satisfies y > alogcosh(t/a) in (12. 4p . This 
means that this string pair is annihilated for y <^ alogcosh(t/a). They are recombined 
at y = alogcosh(t/a). In summary, this D-brane at a fixed time looks like a hairpiiu as 
described in the Fig{Tl 



^In the presence of a linear dilaton, the hairpin shaped D-branc is ahowcd as a static configuration 
^^ as opposed to our case. 
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Figure 1: The left figure shows the hairpin-like shape of the time-dependent D-string 
(12.41) at a fixed time t. The right one describes its time-dependence explicitly. The tip 
(or turning point) of the D-string is located at x = and its time evolution is given by 
y(t) = alogcosh(t/a). This trajectory becomes light-like in the limit t -^ ±oo. 



The effective action of curved D-branes is described by the DBI action (12. ip only if we 
can neglect the higher derivative terms. This requires that the curvature of the induced 
metric is smaller than the string scale. As we will briefly discuss in the next section, 
such stringy corrections become important only if we consider the behavior near t -^ ±oo 
in the region close to the recombined point (or turning point) x = 0. Thus the stringy 
corrections are not essential to our main results in this paper. 

One may notice that another possible decay process is that the brane and the antibrane 
approach via the attractive force between them. A simple analysis estimates the life time 
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due to this process as ~ "3 . This is clearly much longer than the time scale a of our 

•/SSI's 

solution ([23D. 

We can take into account the gauge field A^ on the brane. A solution with an electric 
flux is shown in the in the appendix El We can also generalize straightforwardly the above 
analysis to the systems of a Dp-brane and an anti Dp-brane by trivially adding extra 
coordinates. The profile (12. 4p remains the same in this case. In these higher dimensional 
examples, there are also other decay modes such as making a hole leading to a wormhole 
like configuration as discussed in [12] . We leave the analysis of the time evolution of such 
configurations for an interesting future problem. 



3 Induced Metric and Transverse Scalars 

Open strings on the D-string probe the time dependence of the background via the 
time-dependent induced metric of the D-string 

ds^ = „ H — — -— + 2 tan(x/a) tanh.(t/a)dtdx. (3.1) 

cosh (t/a) cos^(x/a) 

3.1 Induced Metric in Conformal Gauge 

It is useful to rewrite this metric into the one in the conformal gauge 

ds^ = A{T,y){-dr^ + dy^), (3.2) 

so that the Lagrangian of minimally coupled massless scalar fields becomes simple at 
quadratic order 

L = l{-idr<t^f + idy<Pr)). (3.3) 

Requiring that the coordinate y in (13.21) is the same as the one in (12. 4p . we find an 
appropriate coordinate transformation as follows 



r = a log Y 1 + sinh (t/a) sin^(a;/a) + sin(x/a) sinh(t/a) 
y = alogcosh(t/a) — Q;logcos(x/a). (3.4) 

It is equivalent to the relation 

sinh(r/a) = sinh(t/a) sin(x/a), e^^/" = ^"^^(^Z") . (3.5) 

cos(x/a) 

The metric in this coordinate system reads 

, „ (1 cosh^ \ , 

ds^ = - \- + — ^° = dudv, (3.6) 

y 2Jsinh^-sinh^y 

where we defined 

u = T — y, V = T + y. (3.7) 

We need to deal with the metric (13. 6p carefully. It is well-defined for y > \t\. When 
we go beyond this bound, the coordinate y (and r) in (13.61) becomes time-like (space-like) 
and we have to exchange r with y and vise versa, which is equal to the sign flip of u 
and V. Also at the bound y = ±r, the metric becomes singular. Actually, we can further 
perform a coordinate transformation and obtain a completely smooth metric which covers 
the total spacetime on the D-string as we show in the appendix [Bl 



3.2 Asymptotically Flat Coordinate 

The metric (13.61) at the null past infinity u = —oo approaches 

A(u,oo) = - (l+ , l V (3.8) 

On the other hand in the limit v = oo oi the future null infinity we get 

A(-oo,v) = - (l + ^= I . (3.9) 

Thus we can define the asymptotically fiat metric by the coordinate transformation 



a 



V = log(l - VI - e-V"), 

Ct I 

f/ = -log(l- Vl-e2«/"). (3.10) 

In terms of (f/, V^), the metric ds^ = —A{U, V)dUdV is found to be 

A(U,V) ^ Aiu,.)^^ (3.11) 

2 fl + e-W«e^/"V2 - e-2V/aV2 - e2^/° + (1 - e-2W")(l - e^^/") 



(2-e-2W°)(2-e2f^/°) 
It is easy to see that this is manifestly asymptotically fiat: 

ds^ ~ -dUdV (3.12) 

near the past null infinity {U -^ — oo) and the future null infinity (V -^ oo). This 
coordinate system is regular in the regime 

y>--ln2, U<-\n2. (3.13) 

3.3 Global Geometry 

In this coordinate {U, V), the turning point x = is described by the curve 

e^^+e— =2, (3.14) 



or 

a 



V = P{U), P(t/) = --ln(2-e2^/"). (3.15) 

The half of the D-string spacetime x > is mapped to the region V > P{U), while another 
half a; < is mapped to the region U > P{V). Since they are actually equivalent by the 
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'U=tanl/2(\|/-^) 
.V=tan 1/2(\|/+^) 




V=- a/2 log 2 



U= a/21og 2 V= «/21og 2 



V^ oo U^ oo 



71^ -^ 




U=- a/2 log 2 



Figure 2: The Penrose diagram of our time-dependent spacetime (13.1 II) . a; > is mapped 
to region I, and x < is mapped to region II. The global geometry can be obtained by 
gluing two shaded regions along the curve x = represented by the thick curve. 

reflection, the global geometry is obtained by gluing two copies of the region V > P{U) 
along the symmetry axis V = P{U) as in Figj^l Note that both of them are included 
in the regime (13.131) of validity of the coordinates (f/, V). For more explanations of the 
global structure of our metric, see the appendix [Bl 

If we deflne U = i — y and V = i + y, this curve is described by 



— log cosh — . 
2 ^ a 



(3.16) 



Notice that this curve is slightly off from the naive estimate y ^ a log cosh - from (12. 4p 
by the factor 2. This is related to the fact that the value of the metric function A on 
V = P{U) is off from unity, 

i(f/,P(f/)) = 4. 



On the other hand, along the symmetry axis dU and dV are related as 

f,-2V/a 



dU 



-dV. 



2 - e-2^/" 
Thus, the induced line element along the symmetry axis is 



(3.17) 



(3.18) 



ds' 



-2V/a 



dV\ 



(3.19) 



This shows that the symmetry axis V = P{U) is asymptotically null in the limit V/a -^ 

CXD. 

One can check that the extrinsic curvature of the symmetry axis vanishes so that the 
two copies of the region V > P{U) are glued smoothly. For this purpose, let us introduce 
a unit normal vector n" as 

rf = Le-f^/«e-^/", r/ = ^e^/"e^/". (3.20) 



The extrinsic curvature is defined as 



Kab = T^Qaqifi^cnd + V^nJ 



(3.21) 



V=P{U) 



where 

Qab = gab - UaUb. (3.22) 

By explicit calculation it is shown that 

Kab = 0. (3.23) 

Therefore, two copies of the region V > P{U) are glued smoothly as should be. Hereafter, 
we call one of the two copies of the region V > P{U) the region I and the other the region 
II (see Fig.2.). In this way we get a completely smooth coordinate which covers the total 
spacetime. 

Before we proceed, we would also like to examine stringy corrections which we ne- 
glected in this paper. In order to argue it is neghgible, we need to require /^-R/jjpo- <^ 1 
in the region I (the analysis of region II is the same). The curvature tensor takes larger 
values as we come close to the curve V = P{U) and there it behaves like R\iupa ~ ^e^^/°. 
Therefore the condition of neglecting stringy corrections is given by 

V^«|log^. (3.24) 

Since we assumed from the beginning a ^ Z^, we can still neglect the stringy corrections 
even if we take — to be appropriately large. 

3.4 Transverse Scalar as Tachyon 

The main motivation of the above analysis of the induced metric is to understand the 
low energy effective action of open string modes. Since gauge fields in 1 + 1 dimension 
are not dynamical, we would like to consider the transverse scalar fields. First we focus 
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on the transverse scalars </)j (z = 3,4, ■ ■ -, D) which correspond to the deformations into 
the directions orthogonal to the R^'^ defined by {t,x,y). In the low energy limit, their 
effective actions are all equal to the one of the (minimally coupled) massless scalar field 
-^ = ^V~gg^'^9^(j)di^(f), where /x and u run over the D-string world- volume and g^y is the 
induced metric of the D-string (13. lip . 

There is another transverse scalar which describes the motion of D-string within the 
-R^'^ spacetime. As shown in the appendix [Dl at the quadratic order, the effective action 
is given by the one of a massive scalar field L = ^y/—g{g"'^da(j)db(l) + rri^cj)^) with the mass 
square given by 

m' = ]^{R'^^'^-K>^'K,y), (3.25) 

where K^^ is the extrinsic curvature. It turns out that in our example of the decaying 
D-string and anti D-string the mass square becomes negative 

„ 2sec^- sech^- 

177,2 a a ^ g /g 2g) 

a2(sec2 ^ sech^^ + tan^ ^ tanh^ ^f ' ^ ' ' 

^ a a a a' 

and therefore the scalar is a tachyonic mode. In this paper we restrict to the case where 
the condensation of this tachyon does not occur. It is a quite interesting future problem 
to see what will happen if the tachyon condenses. 



4 Energy Flux from Decaying D-brane 

Now we move on to the main part of this paper. We would like to analyze the particle 
creation (open string creation) in the time-dependent open string theory on the D-string 
(12.41) . As we will mention, this analysis is essentially the same as that of moving mirror 
systems. First we concentrate on massless transverse scalarCl and later we will discuss 
the result when we take all stringy modes into account. 

4.1 Energy Flux in Massless Scalar Theory 

We introduce a minimally-coupled, massless scalar field 0. By decomposing the field into 
left- and right-moving parts as 



h+{V)+h-{U) (region I) 
(t>n+{V)+(t>iUU) (region II) ' 



^We will not discuss the tachyonic scalar (|3.26p because its mass is time-dependent and the quantum 
analysis will be much more complicated. 
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the smoothness of the field imphes that 

h4u) = hi+{Piu)), hi-{u) = h+{P{u)). 

Thus, we obtain a general solution 

h+{V) + hi^PiU)) (region I) 
0,,+ (r) + 0,+ (P(t/)) (region II) 

We consider the vacuum \in) defined by the following mode functions 



1 _Q-ii^V 



/i-KUJ 
_J_ 



(region I) 
^g-i<^p{c/) (region II) 

region I) 
region II) 



1 ^^-iu>p{u) (region I) 



/i-KUJ 



1 .g-JwV 



/i-KUI 



The Weight man function in each region is 

G{x,x') = {in\(f){x)(j){x')\in) 

) 

duj 



£H^)0£^*(^') + 0!fH^)0£^^') 



G^^\x,x') + AG{x,x') 



where 



and 



G(°)(x,x') 



du 

4:71UJ 



^_,^(V-V") _^ g-ia;(C/-C/') 



AG{x,x') = G{x,x')-G'^°\x,x') 

"°° du 







4:TTUJ 






^-iu>{P{U)-P{U')) _ ^-iuj(JJ-U') 



U -U' + ie 



P{U) - P{U') + le 

Here, to obtain the last line, we have used 

f^ duj 
Jo 



UJ 



^gia;x _ ^iu^y^^ ^ ^^ 



y-ie 



X — le 



(4.2) 



(4.3) 



(4.4) 



(4.5) 
(4.6) 



(4.7) 



(4i 



Thus, the renormalized stress-energy tensor for each field for the state \in) is decom- 
posed as 



-i(ren) rTi(0) 



■ ah 



' ah 



AT, 



ab; 



(4.9) 
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where T^^ is the contribution from G^^^ obtained later and AT^b is the contribution from 
AG given by 

ATuu = \imdudu'AG{x,x') = --^{P{U),U}, 

X'^X 24:71 

ATvv = limdvdv'AG{x,x') = 0, 

x'^x 

ATuv = ATvu = 0. (4.10) 

where {/, z} is the Schwarzian derivative defined by 



if,r.}^-2Vm^-l-^=], (4.i; 



and satisfies 



The contribution T^\, from G^^\x,x') can be obtained by properly subtracting di- 
vergent parts, differentiating with respect to x and x' and taking the coincident limit. 
The result must be local in the sense that all components are functions of the metric 
function A and its derivatives. In particular, they must be independent of the function 
P{U). Moreover, T^^ must vanish in the asymptotically fiat region since the subtracted 
divergent parts are local geometric functions. 

Actually, all components of T^^ can be determined by the following three conditions: 
(i) it satisfies the conservation equation V"T^^ = 0; (ii) its trace is determined by the 
trace anomaly t'^^" = i?/(247r), where R = 4A-^[AdudvA - {duA){dvA)] is the Ricci 
scalar; (iii) it vanishes in the asymptotically fiat region. The general solution to the 
conditions (i) and fii) is 



-'-uu 






Tlji^ = TZ = -^R, (4.13) 

where a{U) and /3{V) are arbitrary functions. Finally, the condition (iii) implies that 

a{U) = 0, (3{V) = 0. (4.14) 
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In summary, we obtain 



p{ren) 

■uu 



rp{ren) 

J-vv 



12n 

1 
12^ 



p{ren) 
UV 



where 



n(0) 



AT, 



Adl 



Adl 



2AtiA^ 



AT, 



uu, 



uu 



247r 



{duA){dvA) - AdudvA 



{PiU),U}. 



(4.15) 
(4.16) 



In the asymptotically fiat region, only ATuu is relevant: 



rpiren) 
rp{ren) 

J-vv 

rp{ren) 



AT, 



uu, 



n(0) 

vu 



0. 



For P{U) in (I3.15p . the Schwarzian derivative is 

8 1 - e^^/" 



{P(f/),f/} 



Thus, we obtain 



AT, 



a2(2-e2^/°)2' 



1 1 - e2f^/" 



uu 



(4.17) 



(4.18) 



(4.19) 



37ra2(2-e2f^/")2" 
Notice that this energy flux (I4.19P approaches the following constant in the limit 

ATuu - TTT^ = ^Tl (4.20) 



127ra2 12 ^ 

This suggests that a static detector will observe the thermal flux at the temperature 
Td = — long before it will collides with the tip of the D-string (i.e. x = 0). This can be 



more clearly conflrmed by computing the Bogolubov coefficients as we do in the appendix 

o 

However, later it gradually decreases and eventually changes its sign (i.e. negative 
flux) as described in Figl31 In the end the detector collides with the tip of D-string. Thus 
the negative flux is observed only for a flnite time At = ^ log 2. 

It even gets negatively divergent precisely when it collides with the tip of the D-string 
at t/ = flog 2 if we literally accept (14.191) . However, we have to remember that the 
string corrections become important when the condition (13.241) is violated. Indeed, if 
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U=a/21og2 



U=0 



positive ilux 



-? 

negative flux I 




thermal flux 



Figure 3: The behavior of the energy flux due to the particle creation under time evolution. 
In the left figure we show the spacetime regions where the sign of the energy flux is + 
or — . A detector which is far from the tip of the brane first observes a thermal energy 
flux. However, later it changes into a negative energy flux and finally collides with the 
end point. This is summarized in the right figure. 

we assume V ~ |logp, we find that the energy density becomes of the same order as 
the fundamental string tension ATj/jy ~ p-. Therefore when this divergence happens the 
stringy corrections cannot already be negligible and they will plausibly give substantial 
modifications. 

4.2 D-brane Interpretation and Moving Mirror 

Now we would like to interpret the above result from the viewpoint of the decay of the 
D-string-anti D-string system. Before the decay occurs, the effective world volume theory 
of the system at low energy is regarded as a 1 + 1 dimensional f/(l) x t/(l) gauge theory. 
Concentrating on a massless scalar in this theory, the ones for the brane and antibrane 
are denoted by 0/(t, y) and (pnit, y). Here notice that we considered that the fields on the 
D-brane and the anti D-brane live in the same 1 + 1 dimensional spacetime assuming that 
the distance between the D-branes is much smallecl than the typical length scale which 
the observer probes. 

After the time evolution (12. 4p starts, the brane and anti-brane become connected at 
a point y = y^:{t) and the U{1) x f/(l) gauge symmetry is broken there. The continuity 



^We still assume that the distance is much larger than string scale. 
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of the scalar fields at y = y*{t) requires the following boundary conditions 

Mt,y*{t))-<Pii{t,y*{t)) = 0, n''da[Mt,y) + Mt,y)]\y=y^(^t) = o, (4.2i) 

where n" is a unit normal to the curve y = y*(t). Taking the global metric, the function 
y*{t) in our case is explicitly given by (13.161) . 

In our global geometry, the massless particle which propagates in the x < (or x > 0) 
region "reflects" at y = 2/*(t), and then comes into the a; > (or x < 0) region. Thus 
we can consider the trajectory y = yx{t) as a perfect "mirror" in the 1 + 1 dimensional 
space-time, and this setup is essentially the same as the one of the moving mirror system 
in which boundary conditions depend on time (see e.g.fTIH fTT]). 

In general situations, when the mirror is moving at a non-zero acceleration, the static 
observer detects the thermal flux by the Unruh effect. It is described by a quantum 
field theory in the presence of the time- dependent boundary condition. For example, the 
Hawking radiation from a Schwarzschild black hole can also be regarded as a particu- 
lar setup of moving mirror defined by y{t) = —t — e"^'^"^'''^*, where T^h is the Hawking 
temperature of the black hole. A similar analysis of the energy flux ATuu in section 3.1 
and of the Bogolubov coefficients in the appendix [Cl correctly reproduces the well-known 
properties of the Hawking radiation [10] . 

Indeed, our boundary conditions (14.211) are equivalent to the moving mirror system 
with two massless scalars. One of the two scalars obeys the Dirichlet boundary condition 
and the other does the Neumann one. Therefore, the energy flux in the low energy world 
volume theory of the brane-antibrane system should be identifled with the twice of ( 14.19p . 
If we consider other massless flelds, assuming the conformal invariance, we just need to 
multiply the value of central charge c with (14.201) . 

In the actual open string theory, there are inflnitely many massive flelds in the brane- 
antibrane system in addition to the massless flelds we discussed. Let us try to take all 
stringy massive modes into account in our argument of particle creation. Notice also 
that if we are interested in higher dimensional Dp-branes {p > 1) case, we can regard 
the transverse momenta as the Kaluza-Klein momenta, generating extra mass terms. 
Therefore the treatment of decaying Dp-branes is similar. Unfortunately, even in two 
dimension, an exact analysis of a massive scalar fleld in the moving mirror setup is not 
available at present. Still we can roughly estimate the total energy flux as follows. 

We consider the open string creation long before the observer collides with the tip of 
the brane. The result (14.191) argues that we will get a thermal radiation at the temperature 
To = — . We can show this explicitly by computing the Bogolubov coefficients from the 
reflection condition (I3.16p as in the appendix Oof this paper. The string production rate 
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is roughly estimated by multiplying the Hagedron density of states {Th is the Hagedron 
temperature) 

dE e^e~^. (4.22) 



Notice that open string modes become in general massive, but their masses clearly sup- 
press the particle creation and so the above calculation gives the upper bound. Thus 
we can conclude that the total string creation is (UV) finite since Th ^ j- ^ '^d under 
the assumption. This result makes a sharp contrast with the open and closed string cre- 
ation in the rolling tachyon backgrounds [151 [HI [13 [21 [H] , where the string creation rate 
diverges. 

4.3 Analogue of Unruh Vacuum and Negative Energy Flux 

In the actual decay processes of brane-antibrane systems, the time evolution like (12. 4p 
will be triggered by tunneling processes or by collisions with other objects such as closed 
strings or D-branes. From this viewpoint, it is more natural to assume that when t < 
the D-brane remains static and that the real time evolution (12.41) starts at the time t = 0. 
We may think that this is analogous to the Unruh vacuum in an analogy to the Hawking 
radiation from black holes. We can equally treat this case using the conformal mapping 
technique. The upshot is that the energy fiux for f/ < becomes vanishing and that there 
remains the negative fiux when f/ > as in FigHl This can easily be understood from the 
causal evolution. If we want to explicitly confirm this, it is useful to start with a smooth 
time-dependent evolution which is very similar to this. For example, we can choose the 
trajectory of the turning point as 

^ = ilog(e^ + l). (4.23) 

a Z 

Notice that now our new y converges for t <^ as y ~ 0. We eventually find the following 

energy fiux 

1 1 -2cosh[2M/a] 

12Tia^ sinh(2M/a) 

Notice that this fiux is always negative and Tuu diverges as m ^ as expecteqj. A 
similar negative energy fiux was found in the analysis of c = 1 matrix model dual to two 
dimensional string theory [18]. 

The negative energy fiux may seem counterintuitive at first. In fact, this phenomenon 
is known to be typical in the moving mirror systems and a number of works have been 






^[/ is related to u via w = [/ — ^ log 2. 
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U=a/21og2 
^^ U=0 



zero flux 



Figure 4: The behavior of the energy flux in the vacuum analogous to the Unruh vacuum 
under time evolution. Notice that the energy flux is always non-positive. 



done to check that it does not violate fundamental laws in thermal dynamics and so on 
[TOl [TT] . A more familiar example of negative flux will be the evaporation of black holes, 
which can also be explained as the moving mirror effect. 

In our setup, observers situated far from the tip will lose their energy due to the 
negative energy flux for a while before they collide with the tip. Since the total energy 
should be conserved, this energy lost by the observers will be gained by the moving D- 
string. 



5 Conclusions and Discussions 

In this paper we analyzed an explicit decaying process of the brane-antibrane system 
described by the simple solution (12. 4p to the DBI action. At a fixed time, its world volume 
takes a hairpin shape (Fig{T]). The position of the tip (or the turning point) in this hairpin 
is time-dependent. This will be one of the simplest time-dependent backgrounds in string 
theory. 

After we calculated the induced metric of its world volume, we computed the quantum 
energy flux due to the particle creation in this time- dependent system. We found that 
it looks like a thermal flux when an observer is far from the turning point. At a later 
time, slightly before the observer collides with the turning point, the energy flux becomes 
negative. This means that observers lose their energy before they collide with the tip of 
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the D-brane. It is an interesting future problem to investigate the physical consequences 
of this phenomenon. Since the energy flux computed in this paper comes from the open 
string excitations, it will be intriguing if we can also compute closed string radiations 
from the decaying brane system. 

Even though the solution itself (12. 4p describe the decay of a Dp — Dp for any p, we 
mainly restrict the analysis of energy flux to the simplest case p = 1 for simplicity. The 
study of higher dimensional case p > 1 will also be an interesting future problem even 
from the viewpoint of phenomenological applications. 

In our analysis, we always described the dynamics of D-branes by the DBI action. 
Therefore we neglected the stringy effects. Even though these corrections are not impor- 
tant for the regime we are mostly interested in, they become relevant in order to know 
what will happen just before the observer collides with the tip. Therefore, it would be 
very intriguing to construct corresponding boundary states which are a' exact. 

When we calculated quantum energy flux from inhomogeneous decay of a brane- 
antibrane pair, we treated the solution (12.41) as a fixed background and totally neglected 
backreaction due to the flux. Once the backreaction is taken into account, one should 
then be able to answer various physical questions. One of such questions would be about 
consequences of the late-time negative energy flux. This is, in a sense, similar to ques- 
tioning physical consequences of quantum energy flux from a black hole. It is well known 
that quantum fields in a fixed black hole background exhibit thermal flux called Hawking 
radiation. Taking the backreaction of this positive energy flux into account, it is con- 
cluded that the black hole should lose its energy, i.e. mass, and that the flux from the 
hole should increase since the hole's temperature is inversely proportional to its mass. 
This catastrophic loss of the black hole mass should continue until the system exits the 
regime of validity of the low energy effective field theory. 

In the case of the decaying brane- ant ibrane system, the late-time energy flux is nega- 
tive and this should correspond to increase of background energy. Therefore, it is naturally 
expected that the motion of the decaying brane- ant ibrane pair should receive extra ac- 
celeration as a result of the negative energy flux. (Note that, even without including 
the backreaction effect, the system receives acceleration since it loses the rest mass of 
the decaying segment. The backreaction due to the negative energy flux adds further 
acceleration to the motion of the system.) In particular, if we consider corrections to 
the behavior shown in the right panel of Fig{Tl the tip of the bent D-string should be 
further accelerated towards the +y direction at late time. It should be an important and 
interesting future work to investigate the backreaction of the negative energy flux more 
quantitatively. 
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Finally, from the string theory viewpoint, it is very likely that the beginning of our 
universe includes many D-branes and anti D-branes. This makes the study of their anni- 
hilation processes very important to understand the cosmology of our universe. Therefore 
it will also be a quite interesting future project to consider applications of our setups 
and results to cosmological scenarios. One problem in this direction will be to analyze 
a similar time-dependent D-branes in AdS or warped backgrounds. Indeed, as stated in 
the introduction, the end of warped brane inflation is annihilation of a brane-antibrane 
pair in a warped geometry, and inhomogeneities of the type considered in this paper are 
expected to be crucial in the (p)reheating process after inflation as well as in the evolution 
of cosmic string network. The inhomogeneous and time-dependent background presented 
in this paper and its extensions to higher dimensions and/or to warped backgrounds are, 
therefore, of cosmological interest in the framework of warped brane inflation. Further 
studies in this direction are certainly worthwhile. 
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A Solution with Electric Field 

Let us take into account the gauge field A^ on the brane and construct the exact 
solution for DBI action as we did in section 2. The electric fiux is given by / = Ft^^ = 
dxAt — dfA^. The lagrangian looks like 

L = Vl- 2/2 + (!/')'-/'• (A. 1) 

The equations of motion reads 
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Thus we can set 7; = c, where c is a constant. Since L = .^ ,^ \/\ — y"^ + y''^, the first 
equation in flA.2p is the same as the one without any flux. Thus we find that the solution 
is given by 



t X 

y = a log cosh a log cos — , 

a a 



f 



c 



sinh 



2 t 



sm 



2 X 



1+C2 



cosh^ - cos^ ^ 



(A.3) 



B Details of Induced Metric on the Moving D-brane 

Here we present the details of the analysis of the induced metric and the useful 
coordinate transformations employed in section 3. In particular we construct a smooth 
coordinate {U, V) which globally covers the total spacetime. 

We start with the metric (13. ip and rewrite it into the expression in the conformal 
gauge. Requiring that the coordinate y is same as the one in (12. 4p . this is done by 
performing the following coordinate transformation 



aArctanh 



tan(x/a) tanh(t/a) 



1 — tanh (t/a) + tan^(x/a) 



a log 



1 + sinh^(t/a) sin^(x/a) + sin(x/a) sinh(t/a) 
y = alogcosh(t/a) — alogcos(a;/a). 
Equivalently, it is also expressed as follows 

sinh(r/a) = sinh(t/a) sin(x/a), e^'" 

Note that it is also possible to express t in terms of r by 

=^/" I cosh(y/a) 



cosh(t/a) 
cos(x/a) 



cosh (t/a) 



smh ■ smh 



a 



a 



(B.l) 



(B.2) 



:b.3) 



This coordinate transformation fIB.ip leads to the metric, which is manifestly confor- 
mally fiat 

ds"^ = A{T,y){-dT'^ + dy^), (B.4) 

where the conformal factor is given by 

1 1 cosh{y/a) 



A{T,y) = l + 



tan^(a;/a) — tanh^(t/a) 



2jsinh^±^-sinh^^ 

a a 



:b.5) 
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After we define 

u = T — y, V = T + y, (B.6) 

the metric is expressed as l \3.6\i . 

We need to deal with the metric fl3.6p carefully. It is a good metric for y > \t\. When 
we go beyond this bound, the coordinate y (and r) in (13.61) becomes time-like (space-like) 
and we have to exchange r with y and vise versa, which is equal to the sign flip of u 
and V. Also at the bound y = ±r, the metric becomes singular. Actually, we can further 
perform a coordinate transformation and obtain a completely smooth metric which covers 
the total spacetime on the D-string. 

This can be done by considering the following coordinate transformatioio 

a a 

Here we took the absolute values of u and v and this corresponds to the exchange of time 
and space coordinate as we mentioned. Equally, 

U = -a log (el"!/" + Ve^l"!/" - l) , V = alog (el^l/" + ^/e^W'^ - l) . (B.IO) 

In this global coordinate system, the metric reads 

ds"^ = -A{U, V)dUdV, (B.ll) 

where 

- ~ ~ tanh — tanh — / 1 + cosh — cosh — \ 

A(U, V) = ^ ^ • 1 + '^ " . (B.12) 

^ \ smh — smh — ^ / 

The above metric is regular at m = and v = 0. Thus this metric covers the total space 
time and U and V take the values —oo<U,V<oo. In this coordinate system, the 
turning point a; = is mapped to U = V, while the originally singular null lines uv = 
is mapped to the smooth null lines UV = 0. 



^Another way to reach this coordinate transformation is to change the metric at t = into the one in 
fiat space 

ds^ = -dt^ H — - = -di^ + df, (B.7) 

cos"^ - 

a 

where y = a tan — . Then the corresponding hght-cone coordinate defined by 

U = t-y, V = t + y, (B.8) 

reproduced the above one (|B.9[) . 
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The metric f IB.llI) at the null infinity V" = oo approaches 



2 

On the other hand in the limit U = — cxd we get 



J ~ ~ , 1 - tanh ^ 
A{U, V) = ^. (B.13) 



2 
Thus we can define the asymptotically fiat metric by 



. ~ ~ ^ 1 + tanh ^ 
A{U, V) = ^. (B.14) 



U= loffcosh — . (B.15) 

2 2 ^ a ^ ' 

In the same way we definqj 

V a V 

^ = 77 + 77 log cosh-. (B.17) 

2 2 a 

Finally, this leads to the asymptotically fiat metric (13. lip , which is employed throughout 
in this paper. 

C Exact Bogolubov Coefficients for Moving Mirror 

Consider a perfect refiecting mirror is placed in a (1 + 1) dimensional massless scalar 
field theory. Following our setup (13.151) . we assume that its trajectory is described by 
y = a log cosh 4, where we rescaled the parameter a in (I3.16P by factor two as a = ^. We 
define the following function p{u) {u = t — y, v = t + y) as in the standard argument e.g. 

m 

p(u) = -aln(2-e"/"), (C.l) 

or 

u/a 

p'(u) = -— ^. (C.2) 

We have the following two choices of normalized outgoing modes: 



I 



Vi 



^-tu^pw ^^^3^ 



TTUO 

Zf = -yL=e-^^'\ (C.4) 



^The inverse relation is 



2(7 2U 2V 2V , , 

e"— =2e~— -1, e—=2e—-l. (B.16) 
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where 0™ corresponds to the reflection of a standard incoming wave, while 0°"* to the 
standard outgoing mode. In other words, the vacuum with respect to 0™ (or 0°"*) defines 
the one for in vacuum (or out vacuum). 

Then the Bogolubov coefficient is expressed as follows 



«c.,.' = (C, C'*), P^,^' = (C C'*)- (C.5) 

More explicitly, for example 

«^,^' = -^ f rfy(C5,C'* - 5*CC'*) (C.6) 

J t~ fixed 

= t r^ dn(C9„Cf-9„C€'*)- (C.7) 



oo 



Using the formula like 



_^ uj'l [1 + zona + luj'a) 

" "^ (2 — g«/"y^"-ig*<^''"+«/" = •2»('^+^')a ^(-*- + ^^0^)^(1 + iuj'a) 

_^ ^ ' u)Y{\^iu)a + iu)'a) ' ^ ' ' 

in the end, we find the exact Bogolubov coefficients aa;,a;' and (i^^^i 

^ , _ ^ o»(c.W)a r(l + iu:a)V{\ + ia;^Qi) 

2T\\fujuj' r(l + iu;a + ico''a) 

/o ^ oi(^^..')5 r(l + ^^a)r(l - z^^a) 

27rvco'u;' J- (1 + *^" - lixi'a) 

Their absolute value squares are 

I |2 = a (e2a^(a;W) _ 1) 

'""'"'' 27r(cj + u') (e2"- - l)(e2'5-' - 1) ^^^^^ 



l/^^.^'l - 77-7- -TV/ orvT... TTT^ .-2A,r.yN- (C-H) 



27r(cL; - J) (e2°'^'^ - 1)(1 - e-2"'^'^ 
The total creation of the particle with the energy 10 is given by 



A 



N^ = l ^^Uxr-;^^^^3i' (C.12) 



where A is a divergent constanll^°l. This is because |/3(^,(^'P behaves ~ 2^ e^^^'^-i '^^en 
u' is very large. Therefore the particle creation leads to a thermal spectrum with the 
temperature T = ^K. 



""^"This divergence is also natural from the constant energy flux (j4.20p . If we integrates it for an infinitely 
long time, we encounter the same divergence. 

22 



D Fluctuations on Curved D-brane 

Let us consider an n dimensional timelike hypersurface E in an {n + l)-dimensional 
spacetime. Provided that S is sufficiently smooth, it is always possible to adopt a Gaussian 
normal coordinate system in a neighborhood of S so that the spacetime metric is of the 
form 

ds^ = g^ydx^dx" + dw"^, (D.l) 

where x^ (/i = 0, ■ ■ ■ , n — 1) are intrinsic coordinates on S and w is the geodesic distance 
from S. Note that ^^j, depends not only on x^ but also on w. In this coordinate system 
the hypersurface S is represented as w = and the induced metric on S is g^^u = g^u\w=o- 
With the bulk metric and the above Gaussian normal coordinate system fixed, we 
now consider a slightly perturbed hypersurface S' represented a.s w = 0(x). The induced 
metric on the perturbed hypersurface S' is 

^;.= (^M. + 5.05.0) L=^(.)- (D-2) 

By expanding this with respect to and its derivatives up to second order, we obtain 

V = 9>.u + V' V = 2^M.0 + ^;.0' + d^(f^d,(P + 0(0=^), (D.3) 

where K^^, = ldwg^u\w=o and K'^^ = \dlgf,u\w=o- 

Expanding \^—g' with respect to h^^, and keeping terms up to the second order, we 
obtain 

V^ = v^ 1 1 + ^/i + ^ Q/i^ - Ih^'h,^^ + 0{h') I , (D.4) 



where h = g'^'^hnu and h^^ = g^^g'^'^hp^. By substituting the expression (ID.Sp . we obtain 



V^ =V^{l + K<j) + ^g^-'d.cPd^cj) + ]^m'<f + 0{cf) \ , (D.5) 



where 



m" = g'^'^K' +K'- 2K>^''K^,, (D.6) 



K = g^'^Kp^ and K^" = g'^Pg^'^Kp^. 

Now, corresponding to the (n + l)-decomposition of the bulk metric fID.ip . the Ricci 
scalar at w; = is decomposed as 

^^"+'1=0 = ^^"^ - 2^'''^';. -K' + ^K^"K^,, (D.7) 

where w-'^' is the Ricci scalar of g^^. Thus, rn^ in flD.6P is rewritten as 



m^ 



1 {i?W - R^^^^\^^ + K^- K^^Kp,] . (D.8) 
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If the unperturbed hypersurface S is a mininial surface satisfying K = and if R^'^'^^^ 
then this formula is reduced to 

1 



m- = ^{R^"^-K^-'K^,). (D.9) 



2 _ 

~ 2' 
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